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Diversity of coherences and origin of electronic transitions of supermolecular nanoring 
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Quantum coherence is highly involved in photochemical functioning of complex molecular systems. 
Co-existence and intermixing of electronic and/or vibrational coherences, while never unambiguously 
identified experimentally, has been proposed to be responsible for this phenomenon. Analysis of 
multidimensional spectra of a synthetic belt-shaped molecular six-porphyrin nanoring with an inner 
template clearly shows a great diversity of separable electronic, vibrational and mixed coherences 
and their cooperation shaping the optical response. The results yield clear assignment of electronic 
and vibronic states, estimation of excitation transfer rates, and decoherence times. Theoretical 
considerations prove that the complexity of excitation dynamics and spectral features of the nanoring 
excitation spectrum is due to combined effect of cyclic symmetry, small geometrical deformations, 
and vibronic coupling. 


Electronic energy transfer (EET) in organic polymers 
plays a central role in their versatile applications as elec¬ 
tric charge carrying or light emitting/absorbing moieties, 
allowing facile and cheap production of optoelectronic 
devices[1]. Even though various nanomaterials, includ¬ 
ing quantum dots or carbon nanotubes, have been exten¬ 
sively studied, the interest still persists due to the great 
variability of their optical and electronic properties and 
enormous chemical diversity. 

EET in large supra-molecular systems is necessarily 
related to coherent phenomena, arising from the quan¬ 
tum mechanical nature of the inner microscopic con¬ 
stituents. In the conditions of electronic-vibrational res¬ 
onance, quantum mixing of electronic and nuclear de¬ 
grees of freedom of vibronically coupled molecular sys¬ 
tems were suggested to result in a host of phenomena 
beyond the adiabatic approximation[2, 3]. 

Recent development of ultrafast nonlinear spectro¬ 
scopic techniques, such as the two-dimensional electronic 
spectroscopy (2DES) facilitated direct observation of 
quantum superpositions in the form of quantum beats 
(dynamic coherences) [4-10]. Although 2DES combines 
correlation information between excitation and detec¬ 
tion frequencies and simultaneously ensures high spec¬ 
tral and temporal resolutions, the experimental iden¬ 
tification of coherences and determination of their ori¬ 
gin is a very complex issue. The questionable na¬ 
ture of the long-lived dynamic coherences in the iconic 


* These authors equally contributed to this work 

1 These authors equally contributed to this work ; Current ad¬ 
dress: Faculty of Mathematics and Physics, Charles University 
in Prague, Ke Karlovu 3, 121 16 Prague, Czech Republic 

^ Current address: Center for Physical Sciences and Technology, 
Savanoriu Avenue 231, 02300 Vilnius, Lithuania 

§ donatas.zigmantas@chemphys.lu.se 


Eenna-Matthews-Olson (EMO) light-harvesting com¬ 
plex (whether electronic[4] or vibrational[2, 11, 12]) 
still causes debates. Mechanisms related to electronic- 
vibrational mixing[13, 14], inter-pigment correlations[15], 
non-secular quantum transport [16], and inhomogeneous 
broadening[17] were proposed to be important, allowing 
EMO to maintain the dynamic quantum coherence for as 
long as a few picoseconds. 

The role of electronic-vibrational mixing in molecu¬ 
lar systems is of a particular interest, since it was re¬ 
cently shown that diabatic coupling to coherent vibra¬ 
tional modes might enhance the rate of EET[14, 18, 19] 
and of the charge transfer[9, 10]. However, this implies 
existence of both electronic and vibrational coherences in 
the same system at the same time, what has never been 
unambiguously observed. 

Porphyrin nanoring, the chemical structure of which is 
depicted in Eig. lA, was chosen for this quest. It con¬ 
sists of six zinc(H) porphyrin molecules forming a belt 
around an inner hexapyridyl template; the porphyrins 
contain aryl groups (3,5-di(tri-hexylsilane)benzene) at 
meso-positions and are interconnected by acetylenes. 
The fine structure of the peaks identified as Si through 
Sq in its absorption spectrum at 77K (Eig. IB) suggests 
that coherences in 2DES might be clearly resolvable for 
this system. 

Previous quantum chemistry calculations using the 
time-dependent density functional theory estimated the 
lowest-energy Sq —Si transition to be around 10566 cm“^. 
The sequence of almost equally-spaced strong peaks 
(for the nanoring studied here found at 11655 cm“^, 
12295cm“^ and 12862cm“^) were discussed to stem 
from the Eranck-Condon progression of the vibrational 
^ 605 cm“^ mode [20]. 

Two separate 2DES measurements of the fully con¬ 
jugated porphyrin hexamer nanoring were performed at 
77 K using laser pulses, the spectrum of which was cen- 
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Figure 1. Structure (A) and absorption spectrum (B) 
of the porphyrin nanoring. Deformations used in the exci- 
tonic calculation are shown on the right. (B) Measured (black 
solid line) and simulated (red and blue dashed lines) absorp¬ 
tion spectrum at 77 K. Energies of the electronic transitions 
to states Si — Se are indicated by solid vertical lines; vibronic 
transitions are indicated by vertical dashed lines. Frequencies 
of resolved coherent beatings are shown as colored segments, 
connecting the states involved in the corresponding quantum 
superpositions. 
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Figure 2. The absorptive two-dimensional spectra at 
different excitation conditions and at 30 fs of the wait¬ 
ing time. Spectra were obtained using laser pulses centered 
at 880 nm (A) and at 800 nm (B), respectively. In the up¬ 
per panels, laser pulse and nanoring absorption spectra are 
shown. In (A) the signal value is multiplied by a factor of 20 
for the plot range where ui < 11157 cm“^. (C) Zoomed in re¬ 
gion of a degenerate peak in the 800 nm measurement (shown 
by red squares in (A) and (B)) and oscillatory dynamics of 
the “P 34 ” and “P 43 ” peaks. Spectra are drawn using linear 
colour scale, normalized to the maximum of each spectrum. 
Dashed vertical and horizontal lines indicate the energies of 
states Si — Se. 


tered either at 800 nm or at 880 nm, thus covering differ¬ 
ent parts of the absorption spectrum (refer to the upper 
panels in Fig. 2A-B for the corresponding laser spectra). 
Refer to the Supplementary Materials and Methods sec¬ 
tion. 

The 2D spectrum obtained using laser pulses at 880 nm 
covers the lowest electronic transition (Fig. 2 A). It is 
dominated by a strong diagonal peak at ^ 11655 cm“^. 
Other features on the diagonal, related to the absorption 
peaks at 10941 cm“^ (Si) and 11373 cm“^ (S 2 ), are much 
weaker. However, the peaks above the diagonal (“P 12 ”, 
“P 13 ” and “P 23 ”) connecting the three diagonal peaks can 
be clearly resolved in Fig. 2A. Excited state absorption 
shows up as strong negative features below the diagonal 
overlapping with positive peaks. 

Laser pulses with the spectrum centered at 800 nm 
were used to investigate the spectral range of the other 
three most prominent transitions. The corresponding 
2D spectrum is shown in Fig. 2B. The spectrum is very 
rich in features and at least 17 peaks can be clearly re¬ 
solved. Interestingly, a peak on the diagonal at around 


11655 cm“^ consists of two previously not resolved [20, 21] 
contributions separated by ^ 80cm“^. It could be esti¬ 
mated from the position of the off-diagonal peaks, indi¬ 
cated as P 34 and P 43 in Fig. 2C that the energies cor¬ 
responding to these states are around 11600 cm“^ (S 3 ) 
and 11680 cm“^ (S 4 ). 

To determine the energy dynamics via the manifold 
of excited states, the 2DES data were taken as a func¬ 
tion of the waiting time ^2 , thus, making the complete 
3D spectrum. Feature-rich oscillatory evolution has been 
observed throughout the whole 2D area as a 

function of ^ 2 - Decaying dynamics was extracted us¬ 
ing three-exponential-decay fitting with one variable time 
constant of 156-250 fs and remaining lifetimes of 235 ps 
and ::^1 ns. The shortest timescale is related to the down¬ 
ward energy relaxation in the exciton manifold. This pro¬ 
cess is observed as the decay of the peaks on the diagonal 
of the spectrum and simultaneous increase of the peaks 
away from the diagonal if no other competing channel 
exists [22, 23]. The longer timescales represent the relax¬ 
ation from the lowest state of the exciton manifold to the 
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Table L Classification of the observed coherences 
(<700 cm“^). Average of four separate measurements are 
considered and the extracted dephasing times and standard 
deviations (Tt are pointed out. The dephasing times of a few 
coherences could not be extracted with a reliable experimen¬ 
tal error due to the limited resolution. Dephasing times were 
obtained by fitting complex oscillatory signals in the time do¬ 
main (for details, see the Supplementary Text online). 


(jj2 (cm ^) 

Dephasing time (fs) 

Gr (fs) 

Origin of coherence 

±80 

> 500 

- 

Electronic S 3 )(S 4 

±180 

- 

- 

Vibrational 

±235 

360 

160 

Mixed |S:)(S5| 

±307 

280 

130 

Electronic S 2 )(S 4 

±380 

> 600 

- 

Vibrational 

±450 

- 

- 

Vibrational 

±570 

80 

20 

Electronic S 5 )(S 6 

±615 

200 

30 

Electronic S 4 )(S 5 

±695 

190 

40 

Electronic S 3 )(S 5 


ground or the other (for example, triplet) state (for more 
details see the Supplementary Text online). 

Coherent beatings in the 2 D spectra can be conve¬ 
niently visualised by applying Fourier transforms t 2 0 J 2 
to the residuals [24, 25]. The resulting Fourier ampli¬ 
tude dependency is a three-dimensional (cji, co’ 2 , cjs) bulk 
spectrum, the 2D slices (co’i,co’ 3 ), or the so-called oscilla¬ 
tion maps, of which can be plotted for a fixed frequency 
UJ 2 - Notice that evolution of Hermitian conjugate coher¬ 
ences \a){b\ and \b){a\ appears at positive and negative 
co ’2 frequencies, respectively[26]. 

The 2D slices at a few selected 002 frequencies are shown 
in Fig. 3A-D. Judging by their pattern, three types of 
peak configurations can be distinguished, (i) The maps 
at CJ 2 = ±80, ±307, ±570, ±615, ±695, and ±1150 cm-^ 
shown in Fig. 3A-B are diagonally symmetric^ i. e. pos¬ 
itive and negative 002 features are mirror images of each 
other with respect to the diagonal as specifically showed 
in Fig. 3A for co ’2 = ±80cm“^. (ii) Oscillation maps 
at Co ’2 = +380 cm“^ and —380cm“^, shown in Fig. 3C, 
are highly diagonally asymmetric with the features be¬ 
low the diagonal in the co ’2 = —380cm“^ map signifi¬ 
cantly stronger than in the Co ’2 = ±380 cm“^ map. (hi) 
Features in the co ’2 = —235cm“^ and co ’2 = +235 cm“^ 
maps (Fig. 3D) are diagonally asymmetric^ but their am¬ 
plitudes are of the similar magnitude for positive and 
negative frequency maps. Refer to the Supplementary 
Fig. SI for oscillation maps at all other distinguishable 
frequencies. 

We have previously shown that electronic, vibrational, 
and mixed coherences are, in theory, manifested by 
their characteristic patterns and symmetries in the co¬ 
herence maps [3]. The symmetry of the experimental 
maps at co ’2 = ±80, ±307, ±570, ±615, ±695 cm“^, and 
±1150 cm“^ (Fig- 3A and B) indicate that the underly¬ 
ing coherences are of the electronic origin. Dephasing 
times of these coherences are shorter than 300 fs (Ta¬ 
ble I), oscillations are present only in the off-diagonal 
regions in the rephasing 2D spectrum (and only in the 


diagonal peaks in the non-rephasing 2D spectrum; see 
Supplementary Fig. S2), and the peaks at positive and 
negative co ’2 frequencies are of the same amplitude. 

The electronic coherence with co ’2 = +80cm“^ (and 
Co ’2 = —80cm“^) can be assigned to the coherent super¬ 
position |S 4 )(S 3 | (and |S 3 )(S 4 |) of the closely-positioned 
electronic states in the vicinity of 11655 cm“^. The elec¬ 
tronic coherence at co ’2 = ±307 cm“^ shows up only in 
the measurement using laser pulses at 880 nm, signifying 
electronic quantum beats between states IS 2 ) and IS 4 ), 
that could not be excited by laser pulses at 800 nm. 

Beatings with the ±570, ±615, ±695 cm“^ and 
±1150 cm“^ frequencies represent quantum coherences 
|S 5 )(S 6 |, |S 4 )(S 5 |, |S 3 )(S 5 | and |S 4 )(S 6 | (and their Her¬ 
mitian conjugates). We find crosstalk in oscillation maps 
at Co ’2 = ±570, ±615, and ±695 cm“^ due to the over¬ 
lapping of different beating frequencies and limited fre¬ 
quency resolution resulting from the damped dynamic 
coherences. For example, two weaker peaks in the map 
of Co ’2 = ±615 cm“^ indicate the crosstalk from coher¬ 
ences |S 4 )(S 5 | and |S 5 )(S 4 | with ±570 cm“^ frequency 
(Fig. 3B). The crosstalk between different frequencies 
is supported by observation that the weaker peaks are 
shifted away from the co ’3 = co’i ± 615 cm“^ dashed lines, 
parallel to the diagonal. Thus, our findings imply that 
states I Si) through jSe) are of electronic origin in contrast 
to previous assignment [ 20 , 21 ] (it should be noted that 
the previously suggested vibrational progression does not 
follow a displaced harmonic oscillator model). 

We assign beatings with the co ’2 = ±380 cm“^ fre¬ 
quency to the vibrational coherence. This follows from 
the oscillation maps (Fig. 3C), which have the pattern of 
oscillating peaks typical for the vibrational coherence[27]: 
the Co ’2 = —380cm“^ map contains many features below 
the diagonal, while the co ’2 = ±380 cm“^ map is similar to 
the electronic coherence maps presented in Fig. 3A and 
B. In contrast to the electronic coherences, the amplitude 
of the Co ’2 = —380cm“^ map is significantly stronger that 
that of Co ’2 = +380 cm“^ [2, 28]. Detailed analysis of the 
map implies that the vibrational ground state coherences 
|g)(g*| (|g*) denotes some vibrationally hot ground state) 
appear exclusively at co ’2 < 0. Their strong amplitudes 
are therefore related to the long lifetimes of the ground 
state vibrations compared with the ones at electronic ex¬ 
cited states and mapped onto co ’2 = +380 cm“^. Similar 
maps at co ’2 = 180 cm“^, 450 cm“^, and 835 cm“^ (see 
the Supplementary Text online) imply their vibrational 
origin as well. 

However, beatings at co ’2 = ±235 cm“^ point out to 
the mixed coherence, signifying the superposition state 
of the IS 5 ) electronic state and vibronically hot state 
|S|) (380cm“^ vibrational mode). The corresponding 
oscillation map is not typical of neither vibrational nor 
electronic coherences (see Fig. 3D) and the beating fre¬ 
quency 235 cm“^ is equal to the difference between the 
corresponding states. 

To support the assignment of the excited states and 
the natures of the corresponding quantum coherences. 
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Figure 3. Diversity of coherences revealed by oscillation maps with different feature patterns. (A) experimental 
coherence maps of electronic coherences at ijJ 2 = ±80, ±307, ±570, ±695, and ±1150cm“^; (B) experimental and simulated 
oscillation maps for electronic coherence at UJ 2 — ±615 cm“^; (C) vibrational coherence at (jU 2 = ±380 cm“^; (D) mixed vihronic 
coherence at UJ 2 = ±235 cm“^. Intensity of oscillations at each point of a map is indicated by the colour scale; background 
contours show 2D rephasing spectrum at waiting time 30 fs. Each experimental and simulated map for particular frequency is 
independently normalized to the maximal amplitude of either positive or negative frequency map. Dashed lines parallel to the 
diagonal are separated by the value of UJ 2 - Dashed vertical and horizontal lines indicate the energies of transitions to electronic 
states Si through Se- 


two different theoretical models were considered: the 
electronic-only model of six excitonically coupled por¬ 
phyrin molecules and the vibronically coupled dimer 
model. The electronic-only model is used to capture the 
connection between the cyclic symmetry of the system 
and its optical response. The vibronic model supports 
assignment of the vibrational, electronic, and mixed fea¬ 
tures. 

In the electronic-only model, the nanoring is repre¬ 
sented by six two-level sites, representing porphyrins in 
a circular arrangement, each site characterized by a sin¬ 
gle transition dipole, rotated by an angle a = 49° with 
respect to the tangent of the nanoring backbone. Two 
types of deformations—twisting and bending of the ring 
backbone—were taken into consideration (refer to the 
inset in Fig. lA). Due to the full 7 r-conjugation of the 


nanoring, the interaction between any two porphyrins 
cannot be described by the dipole-dipole approxima¬ 
tion and appropriate scaling factors of the coupling con¬ 
stants calculated in the dipole-dipole approximation for 
the nearest, next-nearest, and next-next-nearest neigh¬ 
bors were obtained by fitting the electronic-only model 
to the experimental absorption spectrum (see the Supple¬ 
mentary Text online for more details). The simulated ab¬ 
sorption spectrum of this model is presented in Fig. IB by 
the red dashed line. The agreement with the experimen¬ 
tal absorption spectrum is quite good except of the inten¬ 
sity of the S 5 transition. The absence of vibronic coupling 
in the electronic-only model may be responsible for this 
discrepancy. Particularly, the electronic model provides 
the basis of the observed electronic coherences described 
above. However, additional vibrational/vibronic ingredi- 
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ents are necessary to explain the remaining coherences. 

For qualitative simulations of the vibronic coupling 
in the porphyrin nanoring we used the theoretical ap¬ 
proach developed for a vibronic molecular dimer[3, 29], 
extended by including two vibrational modes. We as¬ 
sume that two electronic states with the highest oscilla¬ 
tor strengths, S 4 and S 5 , are coupled to two vibrational 
modes of 380 cm“^ and 450 cm“^ with the Huang-Rhys 
factors of 0.03 and 0.01, respectively. This allows to sig¬ 
nificantly improve the description of the absorption spec¬ 
trum in the range of S 4 and S 5 peaks (Fig. IB). We also 
performed simulations of the 3D spectra and extracted 
its 2D slices at various 002 values, corresponding to elec¬ 
tronic, vibrational and mixed coherences (Fig. 3B, C and 
D). Although we included only two electronic states in 
the vibronic model, calculated and experimental maps of 
the vibrational coherences at ±380 cm“^, the electronic 
coherences at ±615 cm“^, and the mixed coherences at 
±235 cm“^ are in a very good agreement with the exper¬ 
imental ones, confirming our assignments. 

By this analysis we prove the coexistence of electronic, 
vibrational and mixed coherences in a single system of 
the porphyrin hexamer. Analysis of quantum coherences 
aids in disentangling the energy level structure of the 


excited states and their cooperativity. This turns out 
to manifest even in the absorption spectrum, where we 
identify electronic transitions, the intensity and positions 
of which are non-trivially defined by vibronic coupling 
(together with small ring deformations). In the 3D mea¬ 
surement the whole “zoo” of coherences gets raised and 
they inter-operate to maintain long coherence lifetimes 
and coherent excitation evolution. Such coherent quan¬ 
tum properties of a supermolecular system are reported 
for the first time, but should be general for molecular 
aggregates. 
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